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A new evaluation of the charged Higgs boson decay H + —* W +, y is presented in the context of the 
general two-Higgs doublet model. A nonlinear _R^-gauge which considerably simplifies the calculation 
is introduced and simple expressions are obtained for the fermionic and bosonic contributions. The 
H + — > W +r y branching ratio is analyzed for several values of the parameters of the model. Although 

■^j- ' this decay can have a branching fraction as large as 10 -4 in a certain region of the parameter space, 

it is found that such a region is disfavored by the most recent constraints on b — > 57, g — 2 of the 
muon, Z —* bb, and the p parameter, along with the exclusions from direct searches at the CERN 

£SJ ' e~e + LEP collider. The possibility of detecting this decay at future colliders is discussed. 
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t^- ! I. INTRODUCTION 

<N 

One of the most popular extensions of the standard model (SM) is the two-Higgs doublet model (THDM), mainly 
' because of its simplicity and consistency with the minimal realization of supersymmetry (SUSY) yLj. The THDM 
predicts the existence of five Higgs bosons: two neutral CP-even scalar bosons h° and H a , one neutral CP-odd scalar 
boson j4°, and a pair of charged scalar bosons If The detection of the charged Higgs boson would be an 

irrefutable signal of an extended Higgs sector, and its experimental study might help us to elucidate the structure of 
the Higgs potential and shed light on the group representation of the Higgs fields through the study of its couplings 
to the gauge boson fields. In models which only comprise singlets and doublets of Higgs bosons, the decays of the 
charged scalar boson into a gaug e boson pair, namely H+ -> W + <y [lUfl and H+ W+Z 

are absent at the 

tree level, though they may occur at higher orders. While the decay mode H + — » is forbidden at the tree level 

due to electromagnetic gauge invariance, the H + — > W + Z decay can be induced at this order in models including 
Higgs triplets or more complicated representations Q. In spite of their suppressed branching ratios, these decay 
Oh, modes are very interesting, which is due to the fact that their experimental study may provide important information 
concerning the underlying structure of the gauge and scalar sectors. Apart from being rather sensitive to new physics 
effects 0, these channels have a clear signature and might be at the reach of future particle colliders. 
■ In this work we present an evaluation of the one-loop induced decay H + — > W + j in the framework of the THDM. 
This decay was already studied in the context of the minimal supersymmetric standard model (MSSM) , which is a 
restricted case of the general THDM. At the one-loop, the main contributions to this decay arise from the third- 
generation quarks t and 6, whereas the bosonic sector contributes through the CP-even Higgs bosons h° and H° 
along with the W boson. The fermionic contribution and the effects of a fourth family were studied in Ref. . The 
bosonic contribution was studied in the linear i?^-gauge |3|, which is not suitable for such a calculation since gives 
rise to several complications that can be evaded by the use of a nonlinear i?^-gauge [8|. For instance, while in the 
linear i?^-gauge one has to deal with about 100 Feynman diagrams when computing the bosonic contribution 3], in a 
nonlinear i?^-gauge one only needs to compute 21 diagrams. As for the unitary gauge, the calculation happens to be 
extremely hard due to the terms arising from the longitudinal part of the gauge boson propagators Q. Our aim in 
this work is to reevaluate the H + — ► IU + 7 decay in a nonlinear i?j-gauge, which leads to considerably simplifications 
due to the fact that some unphysical vertices are removed from the interaction Lagrangian. We will see that such a 
gauge not only reduces considerably the number of Feynman diagrams but also renders manifestly gauge-invariant 
and ultraviolet-finite amplitudes. Apart from emphasizing the advantages of using the nonlinear i?^-gauge, we will 
analyze the H + — > W + j decay in some scenarios which are still consistent with the most recent bounds obtained from 
electroweak precision measurements. Indeed, it has been argued quite recently [9j that the parameter space of the 
type-II THDM has become tightly constrained according to the last reported value of the muon anomalous magnetic 
moment a M so we concentrate essentially on the still-allowed region. 

The rest of the paper is organized as follows. A brief review of the model and the nonlinear i?£-gauge used in our 
calculation are presented in Sec. [H] Sec. IIIII is devoted to discuss the most interesting details of the calculation with 
emphasis on the facilities brought about by the nonlinear i?^-gauge. The discussion and conclusions are presented in 
Sec. II VI and Ivl respectively. Finally, the Feynman rules necessary for the calculation are shown in Appendix lAl 
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II. THE CP-CONSERVING HIGGS POTENTIAL IN THE TWO-HIGGS DOUBLET MODEL 

The scalar sector of the CP-conserving THDM consists of two scalar doublets with hypercharge +1: = cj>\*) 
and $3, = {4>2i 02*)- The most general gauge invariant potential can be written as 0] 

V($i,$ 2 ) - /i?($ t 1 $ 2 )+/i2(*2*2)- (/4($i$2) +H.c.) +A 1 ($ t 1 $ 1 ) 2 + A 2 (^$ 2 ) 2 + A 3 ($ t 1 $i)($t$ 2 ) 

+ A 4 ($ t 1 $ 2 )(^$i) + \ (a 5 ($1$2) 2 + (A 6 ($t$i) + A 7 (^$ 2 )) ($1*2) + H.c.) (1) 

It has been customary to impose the discrete symmetry $i — > $i and $2 — > —$2 in order to avoid dangerous flavor 
changing neutral current (FCNC) effects. This symmetry is strongly violated by the terms associated with Xq, and 
A7, but it is softly violated by that associated with ^tf 2 - It is worth mentioning that the last term plays an important 
role in SUSY models. All of the terms in Eq. (JIJ are essential to obtain the decoupling limit of the model, in which 
only one CP-even scalar boson is light. In this work we will assume that this discrete symmetry is softly violated, 
which means that we will take A6 = A7 = throughout the rest of the presentation. 

The scalar potential has to be diagonalized to yield the mass-eigenstates fields. The charged components of the 
doublets lead to the physical charged Higgs boson and the pseudo-Goldstone boson associated with the W gauge field 
as follows 



= -<f>f sin P + <f>f cos P, (2) 
G± = 0±cos/3 + ^sinj8, (3) 

with tan/3 = v%/vi, being V\ (v%) the vacuum expectation value (VEV) associated with $1 ($2)- On the other hand, 
the imaginary part of the neutral components <p®j define the neutral CP-odd scalar and the pseudo-Goldstone boson 
associated with the Z gauge boson. The corresponding rotation is given by 

A" = -4%; sin /3 + <t%j coa /3, (4) 
G z = <t>° u cos p + cfrlj sin p, (5) 

whereas the real part of the neutral components (f^ R define the CP-even Higgs bosons h° and H°: 



h = —(j)i R sina + (j)2 R cosa, (6) 
H {) = 0i_rCoso: + (f)2 R sina, (7) 

where the mixing angle a is given by 



2 7TT.1 2 , \ 

tan 2a = — , (8) 

Tin — m 2 2 

with the elements of the mass matrix being ran = 2vf\\ + tan/3/i^ 2 , Ti22 = 2u|A2 + cot /3/if 2 , and m\i — wi«2(A3 + 
A4 + A 5 ) — fi 12 . 

The masses of the charged and CP-odd scalar bosons satisfy the following expression 



m 2 H± =m 2 A0 + ^f(\ 5 -\ i ), (9) 

where g is the coupling constant associated with the SUl{2) group. Obviously m H ± = m^o when A5 = A4, which is 
a reflect of the underlying custodial symmetry. 

We now would like to comment on the gauge-fixing procedure which was used to simplify our calculation. To this 
end we introduce the following gauge- fixing functions 



./+ = (d; + % -?&Z M ) W+» - OimwG+r , (10) 

f = d^-ZmzGz, (11) 
f A = M". ( 12 ) 
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with D e the electromagnetic covariant derivative and £ the gauge parameter. Note that _f + is nonlinear and transforms 
covariantly under the electromagnetic gauge group. This gauge-fixing procedure is suited to remove the unphysical 
vertices WGwl and WGwZ, which arise in the Higgs kinetic- energy sector, and also modifies the Yang-Mills sector. 
One important result is that the expression for the WW'y vertex satisfies a QED-like Ward identity which turns 
out to be very useful in loop calculations. In particular, in the Feynman-t'Hooft gauge the Feynman rule for the 
Wp (q)W~ (p)A l _ L (k) coupling can be written as 

T^Jf 7 = ~ ie (9^u(p -k + q) p + g vp {q - p)^ + g w {k -p- q) v ) , (13) 

where all the momenta are incoming. It is easy to see that this expression satisfies a QED-like Ward identity. 

The remaining Feynman rules necessary for our calculation do not depend on the gauge fixing procedure. For 
completeness they are presented in Appendix [S] As far as the Yukawa couplings are concerned, we will concentrate 
on the type-II THDM, in which the up-type quarks couple to the first Higgs doublet whereas the remaining fermions 
couple to the second one. 

III. THE H+ -> W + ~i DECAY WIDTH 

In the nonlinear i?^-gauge, the decay H + — » W + -f receives contributions from the Feynman diagrams shown in Fig. 
^ As far as the fermionic sector is concerned, the main contribution comes from the third-generation quarks, which 
induce three diagrams, whereas in the bosonic sector there are contributions from the (H^, <p°) and (W^, <p°) pairs, 
with <f)° = h° or H°. We would like to emphasize that the nonlinear i?^-gauge considerably simplifies the calculation 
of the H + — > W + 'f decay. First of all, the removal of the unphysical vertex W^G^yJ, allows one to get rid of all 
those diagrams which are displayed in Fig. 2 and do have to be calculated in the linear i?^-gauge. Even more, in 
the nonlinear i?^-gauge the tadpole graphs shown in Fig. [3] vanish. Apart from these simplifications, we will be able 
to group the Feynman diagrams into subsets which separately yield a manifestly gauge-invariant amplitude free of 
ultraviolet singularities. 




FIG. 1: Feynman diagrams contributing to the H + — > W 7 decay in the nonlinear iij-gauge. 4> stands for h and H , and 
X+ for H+ and G+ . Compare with Fig. 1 of Ref. 0. 

Once the amplitude for every Feynman diagram shown in Fig. ^was written down, the Passarino-Veltman method 
|12| was applied to express it in terms of scalar integrals, which are suitable for numerical evaluation |13|. The full 
amplitude for the decay H + {p) — > W + {q) r y{k) can be written as 



M(H+ -> W+ 1 ) = — j= T (G(q-k 9lw - kvqj ~iHe^ aP k a qP) e^(k)e v (q), (14) 

which is manifestly gauge invariant. The H function only receives contributions from the quarks: 

H = x l ' Ktb ( 3 SB (H±,w,b.t) + 5{h±,w) (c| C(M) + 2 s % c (t.b))) > (I 5 ) 

°(H±,W) 
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FIG. 2: Feynman diagrams which appear in the linear _R^-gauge but are absent in the nonlinear one. ho, Ho, Ao, H ± , G^ and 
Gz circulate in both tadpoles, whereas W , Z, t and b also circulate in the right-hand tadpole. 




FIG. 3: Tadpole diagrams whose amplitude vanishes in the nonlinear 7?^-gauge. X stands for h°, H°, A , H ± , G%, G z , W ± , 
Z, t, and 6. 



A tb = -^, (16) 

S/3C/3 

where we have introduced the shorthand notation <5( a b j = m 2 a — m\, <5-B( a ,6, c ,d) = B (m 2 , m 2 , m^)—B {m 2 , m 2 , m 2 ^, 1 
and C( ai b) = Co (m 2 H± ,m 2 v ,0,m 2 ,ml,rn 2 ), where B and Co stand for Passarino-Veltman scalar integrals. In addi- 
tion, cp — cos/3, sp = sin/3, etc. From the 8B( abcd -^ definition, it is clear that the H function is ultraviolet finite. As 
for G, it can be written as 

G= ^~Z2 — 72 \Gtb+ { G h±4>° + G G± <jfi + G w±<t>°) J , (17) 

zm H± o (H± w) y ^ =H o, h o J 

where Gab stands for the contribution of the (A, B) pair: 



Gtb = ^tb(m 2 H± (3 (S( t ,b) + C2/3 8(h±,w)) + m w) $B(H±,w,b,t) — 3£(t,&)<5(jf±, w)<5-B(o,i?± 
+ m 2 H ±8(H±,w) (l - (2m? - Cp8( H ±,w)) C(b,t) + (2m? - s 2 fj 8 {H ±, w) ) C (ti6 )) ), 



b,t) 



GR±h° — ^h°H±(8{H±,h )8(H±,w)8B(Q^ WjH ±^ - (2m 2 H± S( H ± h o) + m^ m^) 8Bi yH ±y /H ± h a- ) 

- m 2 H ±8 {H ± tW) (l + 2m 2 H ±C( H ± th0) )^, (19) 



1 When a = 0, m a must be substituted by 0. 
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h° = ( { m whh°, w) + m H± (3m^ - 2m|o)) £B(H±,w,ft°,w) + <5(ff±,/i°) £(h«, wo&E(o,w,h°,w) 

+ m^±5( H ± )M f) (l + 2m% f C{y Vih o)) (20) 

and 

G w ± H o = \w±h ( >(^( m w^(h°,w) + m2 H± - 2m^o - 4m^±)) <5S (/i -± ilv ,ff±,/iO) + 5(H± : w)fi(h° : w)SB( WH ±j l o- ) 

+ m H± S (H±,W) { 2 ^H±,W)- m w) C (ml,,ml ))\ ( 21 ) 

with 



^H±h° = c a-f3 \ 2 S a -j3 m Z H ± + /if 2 + ( S a Sfjtp ) TO^o ] , (22) 

C /3 S /3 



. _ S2(a-j3)^(if±,fe°) 



(23) 



and 



A ff± , = - m ^ (a ^ . (24) 



Finally, the contribution of the heaviest CP-even scalar boson H° is obtained from that of the lightest one once the 
substitutions m^o — * m^o and a — ► a — n/2 are done. 

The above expressions are very simple and should be compared to those presented in Refs. 0,0. It is also evident 
that the partial amplitudes induced by the pairs (t, 6), (0°, W^), (0°, H^), and (0°, G^) are gauge invariant and 
ultraviolet finite on their own. Again, this is to be contrasted with the situation arising in the linear i?^-gauge, where 
showing gauge invariance is somewhat cumbersome, and the cancellation of ultraviolet divergences in the bosonic 
sector is achieved only after adding up all the Feynman diagrams 3] . 

Once Eq. H14[) is squared and the spins of the final particles are summed over, the decay width can be written as 

F(i7+ -> W+,) = 2 J/^ 3 (\G\ 2 + 1*1") • (25) 
z 7r s w m w m H± 

We will evaluate this decay width for some values of the parameters of the model. 



IV. NUMERICAL RESULTS AND DISCUSSION 



In order to evaluate the H + — > W + j decay, there are six free parameters of the THDM which should be addressed. 
They are the masses of the four scalar bosons m H ±, m/jO, m^a and m^o, together with the VEVs ratio tan/3 and the 
mixing angle a. Without loosing generality, we will assume that /ii2 = 0. Although the H + — > W +r y decay width does 
not depend on m^o, the respective branching ratio does depend on all of these parameters. Instead of examining the 
behavior of H + — * W + "f for arbitrary values of these parameters and looking for those which enhance its branching 
ratio, a realistic analysis must take into account the most recent constraints obtained from precision measurements 
and direct searches at particle colliders. We will thus concentrate on the region of the parameter space which is still 
consistent with low-energy data. 

The parameter space of THDMs as constrained by electroweak precision measurements has been recurrently studied 
in the literature 0,^3,01- Along this line, it is well known that the experimental bound on the inclusive decay b — > 
can be translated into a bounded area on the m#±-tan/3 plane. In addition, the muon anomalous magnetic moment 
cifi is useful to constrain the CP-odd scalar boson mass m^o , which along with the CP-even scalar boson masses m^o 
and nifto can be further constrained by the Z — ► 66 decay mode Rb- All these constraints, when combined with 
those obtained from the p parameter and the data from direct searches at the CERN e~e + LEP collider, can yield 



6 



the allowed region in the parameter space of THDMs. Quite recently, the most up-to-date measurements on b — > S7 
E, <V i?b E3. and p together with the DELPHI [H} and OPAL excluded regions, were considered to 
constrain the parameter space of the type-II THDM 0,0]. In particular, the possibility of a very light CP-odd scalar 
boson was considered in |l5j. whereas a more general treatment was presented in From these studies one is lead 
to conclude that the parameter space of the THDM has become tightly constrained by the most recent measurement 
of a M [l3 ■ Below we will follow closely the analysis presented in Ref. [£| • 

The model independent bound m H ± > 80 GeV has been derived at LEP for a charged Higgs boson deca ying solely 
into tv [2]]]. On the other hand, within the THDM-II, the latest reported value for the decay b — > sj |l6j] yields 
the very stringent lower bound m#± > 500 GeV for intermediate and large values of tan j3, though this bound can 
loosen in SUSY models with conserved or broken i?-parity j^. As far as a„ is concerned, the latest experimental 
data from the BNL collaboration [ljj require a positive contribution from new physics to bring the theoretical value 
closer to the experimental one. It turns out that the CP-odd scalar yields a positive contribution to a M , whereas 
the remaining Higgs bosons give a negative contribution. This means that the deviation between the experimental 
and theoretical values of can only be explained by the existence of a ligh t CP-odd scalar along with a high value 
of tan/3 and relatively heavy H^, h° and H°. When the Rb constraint [12] is taken into account, a large region of 
the parameter space of the THDM gets excluded, and only a tiny region survives when the DELPHI 0] and OPAL 
[20| excluded regions for the scalar boson masses are considered. From this analysis, it was concluded that the 
most favorable region is that in which tan/3 is of the order of 50 or larger, m^o < 80 GeV, m h a < 140 GeV, and 
— 1/2 n < a < — 3/8 it. Even more, the p ~ 1 constraint requires fine-tuned solutions for the masses of the scalar 
bosons and the mixing angles. In fact, this constraint is satisfied by an almost decoupling CP-even scalar boson H°, 
with a mass of the order of 1 TeV, provided that the remaining parameters meet the previous constraints. 

Bearing in mind the preceding discussion, we will analyze the H + — > W + j decay for 300 GeV < m H ± < 600 GeV. 
For the sake of illustration we will assume that the lightest CP-even scalar boson is relatively light and the heaviest 
one is much heavier indeed, i.e. we will take m^o = 115 GeV and m^o = 1 TeV. The former value is consistent with 
the current LEP bound on the mass of a SM-like Higgs boson ^^|. Moreover, we will focus on the region where tan/3 
is large and sin a is close to -1. We have verified that these values lie inside the still-allowed region of the parameter 
space of the THDM-II. 
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FIG. 4: H + — > VF + 7 decay width (solid line) as a function of tan/3 for the indicated values of sin a and m H ±. We have taken 
m h o = 115 GeV and m H o = 1 TeV. The partial contributions from the quark (dashed line) and boson (dotted line) sectors are 
also shown. 



It is interesting to examine the effects of both the quark and boson sectors on the decay width T(H + — * W +r y). 
The partial contributions, together with the full contribution, are displayed as a function of tan/3 in Fig. whereas 
in Fig. [5] we show them as a function of m H ±. We can observe that, in this range of parameters, the H + — * W + j 
decay is slightly dominated by the bosonic contribution. The fcrmionic contribution is larger only in the region of the 
parameter space where sin a is very close to -1. This is illustrated in Fig. where we show the partial contributions 
to the H + — > W + j decay as a function of sin a for m#± = 300 GeV and two different values of tan /3. We see that the 
bosonic contribution is maximal when sin a = and minimal when | sina| = 1, where it is smaller than the fermionic 
contribution. For this particular set of the parameters, the maximal value of the bosonic contribution is more than 
one order of magnitude than the respective minimal value. A similar behavior is observed for larger values of m H ± 
and tan/3. Also, from Fig. 0]and Fig. |S]it is interesting to note that the partial contribution from the bosonic sector 
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FIG. 5: The same as Fig. 0]but now as a function of m H ±. 



exceeds the total contribution in some regions of the parameter space. This stems from the fact that the interference 
between the partial contributions is destructive. While the fermion contribution depends only on tan /3 for a given 
m H ± , the bosonic contribution depends on four other unknown parameters. This makes harder the analysis of this 
decay mode. 



T 1.DE-4 




FIG. 6: Partial contributions to the H + — > W +, y decay width as a function of sin a for m H ± = 300 GeV and two different 
values of tan /3. The constant lines correspond to the fermion contribution and the remaining to the bosonic one. The values 
of the CP-even scalar boson masses are the same as in Fig. 0] 

The full width of the charged scalar boson is given to a good approximation by the sum of the partial widths of 
the two-body decay channels induced at the tree-level: H + — ->■ tD, cs, tb, W + A Q , and W + h°. The widths for all these 
decays can be readily obtained and we refrain from presenting the respective expressions here. In the region of the 
parameter space that we are considering, the charged scalar boson would decay mainly into tb. This decay mode is 
somewhat troublesome as it would suffer from large QCD backgrounds at a hadronic collider 23] . The H + — > W + j 
branching ratio is shown in Fig. 0as a function of tan/3 for m^o = 80 GeV and two different values of sin a. We can 
observe that Br(H + — > W +r y) is at most of the order of 6 x 1CP 6 , it increases for larger values of tan/3, but falls slowly 
as sin a gets closer to —1. On the other hand, the behavior of Br(H + — > W + j) as a function of m H ± is depicted in 
Fig. [S] For tan (3 = 40, the branching fraction decreases very slowly for larger mjy± , whereas for tan /3 = 60 it has a 
minimal around to#± = 350 GeV and increases slowly until about m#± = 550 GeV, where there is a maxima. 

It is worthwhile to mention that the H + — > W + ^/ branching ratio may be dramatically enhanced, mainly by the 
fcrmionic contribution, around the threshold m H ± ~ m t + m&. This is shown in Fig. [!|]for some illustrative values 
of tan/3 and sin a. It can be observed that Br(H + — > W + "f) can be larger than 10 -4 when tan/3 is small regardless 
of the value of sin a, and up to 10~ 4 when a large value of tan/3 combines with sin a ~ 0. In the latter case the 
bosonic contribution is also responsible for the enhancement. However, when a large value of tan/3 combines with 
sin a ~ 1, there is large cancellations between the two contributions, which results in a significant suppression of 
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FIG. 7: Branching ratio for the H + —* W + j decay as a function of tan (3 and the indicated values of sin a and m H ± . 
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FIG. 8: Branching ratio for the H + — > W /+ 7 decay as a function of m H ± for the indicated values of sin a and tan/3. 




FIG. 9: The same as Fig. [g]m other region of the parameter space of the THDM-II. 



Br(H + — > W + j). Unfortunately the region of the parameter space in which such an enhancement shows up is not 
favored by the constraints from electroweak precision measurements. Although a fourth-fermion family may improve 
the H + — > W + ^i decay rate ;2], the enhancement would depend on the splitting between the masses of the fourth- 



generation quarks ra u . 



which cannot be too large due to the p ~ 1 restriction. In this respect, it seems that 



the H ± — ► W^Z decay is more promising since it is sensitive to nondecoupling effects when there is a large splitting 
between the CP-odd and the charged scalar boson masses As discussed above, this scenario is still consistent 
with the constraints on the parameter space of the THDM-II. Furthermore, the — * W^Z decay can arise at the 
tree-level in Higgs triplet models . 
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Let us now discuss shortly on the possibility of detecting the H + — » W + pi decay at future colliders. Single and 
pair pr oduction of heavy charged Higgs bosons at both hadronic (23, UM, |25|, |2(| 123, HE H3 an d linear e~e + colliders 
polliulE^ I have long received considerable attention. A charged scalar boson heavier than 180 GeV could be detected 
at a hadronic collider via several processes, such as gb -> H~t |H HI, cs ff+ H3, and c6 -> H+W\. The 
one-loop processes <?<? — * W^H^ and — * W ± H Zf give no viable signature at the CERN LHC collider 

ST While 

the qb — > H~t rate is out of the reach of the Fermilab Tevatron, it could be detected at the LHC for m H ± up to 1 TeV 
[23l |23 ■ This process is appropriate to look for a charged Higgs boson in the low- and high-tan /3 regime. Though 
the gb — > H~t process followed by the decay H~ — > tb suffers from irreducible background, the H~ — > tv mode is 
background free due to the distinctive r polarization [2^. The cross section for gb — > H~t increases for large tan/? 
but falls when m H ± increases. At the LHC, the cross section for pp — ► H + ib can be as large as 10 3 fb for m#± of the 
order of a few hundreds of GeV and tan/3 around 50. With an integrated luminosity of 100 fb, this rate is enough to 
detect a charged scalar boson decaying as H + — > fv, but it is not viable to detect the H + — > mode unless its 

decay width is dramatically enhanced. There is some chance to detect the H + — > W + j decay if the constraints on 
the parameter space can be evaded in some way and the charged scalar boson mass is about 180 GeV. In this case 
the branching ratio might increase up to 10~ 4 , whereas the cross section for pp — > H + tb can reach the 10 4 fb level. 

In a linear collider (LC) the situation is less promising. If the mass of the charged scalar boson happen to be 
smaller than half the center-of-mass energy yfs of the LC, it might be detected through the double pair production 
process e~e + — * H~H + |30|. which has the advantage that it does not depend on tan/3 as it is dictated entirely by 
the interactions of the charged Higgs boson with the neutral gauge bosons. On the other hand, when yfs < m H ± /2, 
a charged scalar boson might be detected via the single production modes e _ e + — > tvH + , e~e + — > bcH + , or 
e~e + — ► W^H^ |23. The first two are more promising as arise at the tree level, whereas the last one occurs at one 
loop. In a linear collider running at y/s = 500 GeV, the cross section for e~e + — > tvH + would be about 10~ 2 fb for 
m H ± = 300 and tan/3 = 60 GeV. Even with an integrated luminosity of 500 fb -1 , this rate is too small to allow the 
detection of the H+ -> IV+7 decay. If y/s = 1000 GeV, cr(e~e+ -> H~H+) - 10 fb for m H ± = 300 GeV, which 
evidently is also too small for the H + — > W + j mode to be detected. 

Since a LC would also work in the mode 77, it has been recently proposed to look for a heavy charged scalar boson, 
with mass larger than half the center-of-mass energy of the LC, through the single production processes 77 — » tvH + 
and 77 — > bcH + J^l). The viability of these modes depends strongly on the strength of the Yukawa couplings of 
the fermions to the charged scalar boson. It turns out that the rates for these processes can be up to two orders of 
magnitude larger than for the respective e~e + processes j^. At a y/s — 500 GeV linear collider, the cross section for 
17(77 ~~ > tvH + ) is of the order of 1 fb, which is a typical value for intermediate rn#± and large tan/3. The detection 
of H + — > W + j would thus not be viable. Of course the charged Higgs boson can also be pair produced in the 77 
mode provided that the respective center-of-mass energy is enough, but the conclusion remains the same. A linear 
collider would thus not be appropriate to look for rare decays of the charged Higgs boson. 

We have mentioned that while the H + — > W + j branching ratio may reach the level of 10 -4 provided that the mass 
of the charged scalar boson is around the threshold m t + nib, it is at most of the order of 10~ 6 in the region of the 
parameter space favored by electroweak precision measurements. Even if the H + — * W + j decay rate was not reduced 
considerably by the kinematical cuts, it would be necessary the production of about 10 6 charged scalar bosons to 
observe a few events. This seems to be beyond the reach of the LHC, so the H + — > W + j decay would hardly be 
detected unless its decay width is dramatically enhanced in other model. 



V. SUMMARY 



A study of the decay H + — > W + j was presented in the context of the type-II THDM. A nonlinear i?£-gauge 
covariant under the electromagnetic gauge group was used to obtain the one-loop amplitudes. The advantages of 
using such a gauge were emphasized. In particular, there are 21 Feynman diagrams in the nonlinear i?^-gauge, 
whereas the linear one induces about one hundred. An analysis of the behavior of the H + — > W + j decay width as 
a function of the parameters of the model was presented, but we focused essentially on the region of the parameter 
space still allowed by the constraints obtained from the latest reported values of b — >■ 57, Rb, and the p parameter, 
which favor a charged scalar boson heavier than the t quark, a large value of tan/3, and —tt/2 < a < —3tt/8. In 
this still-allowed region, which was not considered in the previous analyses of the H + — ► W + 7 decay, the respective 
branching ratio receives contribution of the same order of magnitude from both the bosonic and fermionic sectors, 
and it can be as large as 10 -6 , which seems to be beyond the reach of the LHC. 

Note added: After the submission of this manuscript, the muon g — 2 collaboration announced a new result for the 
anomalous magnetic moment of the negative muon, which is based on the data collected during the year 2001 [33| . 
The new result, a ex f = 11659214(8) (3) x 10~ 10 (0.7 ppm), is consistent with previous measurements of the anomaly 
for the negative and positive muon. The new world average for the muon anomaly is thus a° xp = 11659208(6) x 
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10" 10 (0.5 ppm), which is to be contrasted with the previous one a e * p = 11659203(8) x 10~ 10 (0.7 ppm). This 
new results confirms the fact that THDMs are very disfavored by the muon anomaly. However, we would like to point 
out that the scenario that we are considering is still consistent with the new result. So, the analysis presented in this 
work is not affected significantly and our conclusions remain unchanged. 
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APPENDIX A: FEYNMAN RULES FOR THE DECAY H+ -> W + >y 

In this appendix we present the Feynman rules which do not depend on the gauge fixing procedure described in 
Sec. [n] The H ± H^h Q vertex arise from the Higgs potential. In terms of the coupling constants Xi the Feynman rule 
for this vertex is given by 

9H±K=F/iO = gA°A°h° g^~^ 5 ~ ^) s fJ-a, (Al) 

where 

9A°A«h« = — - — (2 \isj 3 CfjS a - 2\ 2 c 2 pSpc a - (A 3 + A 4 + A 5 ) (s|c Q - cps a ) + 2X 5 Sfj- a ) . (A2) 

The coupling constants Xi can be expressed in turn in terms of the mixing angles and the Higgs boson masses as 
follows: 

2 / / \ 2 / \ 2 

9 I ( s a \ 2 i I C a \ 2 S P 2 



^ = ^{{fe)<° + {fJ m »°-I/»)- (A3) 

^ - 4 ((5) 2 * + fe) 2 '"- (A4) 

2 / 2 
# I o^2 , s 2a / 9 _2 \ Ml2 



A 3 = -V 2mU + —{mU - mfc) - ff- , (A5) 

a 2 / II 2 

y i™ 2 i™ 2 i ^i 2 



g 2 ( \Ai 



^iV -"4 ■ (A7) 

4mf^ \spcp 



After some algebra we are left with 



l g 2 C Q +5 2 / Cj3C a \ 2 1 

gH±HTh» = 7, 2 s a-p m H ± + (4 2 + S a Sfjtp — m h o . (A8) 

imw \ CpSp \ tp ) J 

The Higgs kinetic-energy term induces the H^Gyy h° vertex: 

9 H ±G^ ho = ^;( TO H± - m lo)Cf3-a- (A9) 

The Feynman rules involving the heavy CP-even scalar boson H° can be obtained from above after the replacements 
m/,0 — > m H o and a —> a — it/2. As for the couplings of the fermions to the W gauge boson and the charged scalar 
boson, they are 0: 
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9h- tb 



v + a j 5 ) 



(A10) 



9w~ t b 



2 \f2mw 



(v + a 7 5 ) 7m 



(All) 



with v = m t cot (3 + mi, tan (3 and a = nib tan /3 — m t cot (3. 
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